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Abstract 

The Ruijsenaars-Schneider systems are 'discrete' version of the Calogero-Moser (C- 
M) systems in the sense that the momentum operator p appears in the Hamihonians as 
a polynomial in e^^ ^ (/?' is a deformation parameter) instead of an ordinary polynomial 
in p in the hierarchies of C-M systems. We determine the polynomials describing the 
equilibrium positions of the rational and trigonometric Ruijsenaars-Schneider systems 
based on classical root systems. These are deformation of the classical orthogonal poly- 
nomials, the Hermite, Laguerre and Jacobi polynomials which describe the equilibrium 
positions of the corresponding Calogero and Sutherland systems. The orthogonality 
of the original polynomials is inherited by the deformed ones which satisfy three-term 
recurrence and certain functional equations. The latter reduce to the celebrated second 
order differential equations satisfied by the classical orthogonal polynomials. 



1 Introduction 

Exactly solvable or quasi-exactly solvable multi-particle quantum mechanical systems have 
many remarkable properties. By definition, the entire (or a part of the) spectrum and the 
corresponding eigenfunctions are calculable by algebraic means. The corresponding classical 
systems share also many 'quantum' features. For example, the frequencies of small oscil- 
lations near the classical equilibrium are 'quantized' together with the eigenvalues of the 



associated Lax matrices at the equilibrium. These phenomena have been explored exten- 
sively for multi-particle dynamics based on root systems, in particular, for the Calogero and 
Sutherland systems |2 12 El by Corrigan-Sasaki |4j. Similar phenomena are also reported by 
Ragnisco- Sasaki ^ for Ruijsenaars-Schneider systems [HI 13 IHl E] ? which are deformation of 
C-M systems. 

In this paper we will discuss one special aspect of the classical equilibria of exactly solvable 
systems based on classical root systems, the rational and trigonometric Ruijsenaars-Schneider 
systems. Namely, the determination of the equilibrium positions and their description in 
terms of certain polynomials. It is known that for the Calogero and Sutherland systems, the 
equilibrium positions are described by the zeros of the classical orthogonal polynomials, ie 
the Hermite, Laguerre, Chebyshev, Legendre, Gegenbauer and Jacobi polynomials j^QHlE]- 
The Ruijsenaars-Schneider systems are 'good' deformation of the Calogero and Sutherland 
systems. Here is one interesting evidence. It was known ^21 that the singular vectors of the 
Virasoro and Wn algebras, in the free field representation, are related to Jack polynomials 
[IS], the quantum eigenfunctions of the A type Sutherland systems. The deformed Virasoro 
and Wn algebras were discovered by using the relation between the Sutherland system and 
the trigonometric Ruijsenaars-Schneider system of the A type root system ^Ij. Therefore it is 
expected that equilibrium positions of the rational and trigonometric Ruijsenaars-Schneider 
systems would give certain 'good' deformation of the classical orthogonal polynomials: 
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In Ragnisco-Sasaki paper ^5J, based on numerical analysis, the explicit forms of the lower 
degree members of the one-parameter deformation of the Hermite and Laguerre polynomi- 
als were presented. The present authors continued the numerical analysis and obtained the 
explicit forms of the lower degree members of the one and/or two-parameter deformation 
of the Hermite polynomial, one, two and/or three-parameter deformation of the Laguerre 
polynomial, one-parameter deformation of the Jacobi (and Gegenbauer and Legendre) poly- 
nomials. They are also polynomials, or rational functions in the deformation parameter(s) 
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with integer coefficients. 

Remarkably the orthogonahty of the original polynomials is inherited by the deformed 
ones. The equations determining the equilibrium can be reformulated as functional equations 
determining the polynomials. These functional equations are difference analogs of the cele- 
brated second order differential equations satisfied by the classical orthogonal polynomials. 
Three term recurrence for the deformed polynomials, the necessary and sufficient condition 
for orthogonality, can be derived from these functional equations. Dynamical stability of the 
Hamiltonian system, or bounded-from-belowness of its potentials, is achieved by restricting 
the parameter space of the coupling constants, usually by positive coupling constants, which 
in turn guarantees the positive definiteness of the inner product governing the orthogonal 
polynomials, the deformed as well as undeformed. These deformed polynomials are not the 
so-called g-deformed versions of the above classical polynomials [T3] . 

This paper is organized as follows. In section 2, first we recall the essence of the Calogero- 
Sutherland systems and their equilibria, which are described by the Hermite, Laguerre and 
Jacobi polynomials. Next the Hamiltonians and potentials of the Ruijsenaars-Schneider (R- 
S) systems are briefly recapitulated, and two types of the rational systems and one trigono- 
metric systems for the classical root systems are introduced. Then the equations for their 
equilibrium positions are written down. For later use we review the relation between the or- 
thogonal polynomials and the three-term recurrence. Section 3, 4 and 5 give the main results 
of this paper. In section 3 we determine the equilibrium positions of the rational R-S systems 
for the A type root system and the deformation of the Hermite polynomial is presented. For 
one parameter deformation case, we derive the explicit forms of the generating function and 
the weight function of the inner product. In section 4 equilibrium positions of the rational 
R-S systems for B,C,D, BC type root system are determined and the deformation of the 
Laguerre polynomial is presented together with the explicit forms of the functional equations 
and three-term recurrence. In section 5 equilibrium positions of the trigonometric R-S sys- 
tems for the B,C,D, BC type root system are determined and the deformation of the Jacobi 
(and Gegenbauer) polynomial is presented together with the explicit forms of the functional 
equations and three-term recurrence. Classical orthogonal polynomials, eg, the Hermite 
and Laguerre, satisfy many interesting identities among themselves. Those having a root 
theoretic explanation (folding) are shown to be preserved after deformation. Final section 
is devoted to summary and comments. Identification of the deformed orthogonal polyno- 
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mials within the so-called Askey-scheme of hypergeometric orthogonal polynomials ^E] is 
reported here. The relation between the functional equation and the three-term recurrence 
is discussed in the Appendix. 



2 Potentials and Equilibrium Positions 

In this section we set up models and present the equations for their equilibrium positions. 
We consider a multi-particle classical mechanics governed by a classical Hamiltonian H{p, q). 
The dynamical variables are the coordinates {qj | j = 1, . . . , r} and their canonically conjugate 
momenta {pj \ j = 1, . . . , r}. These will be denoted by vectors in R'' 

9 = *(gi,---,gr), P = *(Pi,---,Pry, 



in which r is the number of particles (and it is also the rank of the underlying root system 

A except for the A case). The canonical equations of motion are 

. _ dH{p, q) . _ dH{p,q) 
^■^ dpj ' dqj 

The equilibrium positions are the stationary solution 

p = 0, q = q, (2.2) 

in which q satisfies 



dH{0,q) 



dq 



0, (j = l,...,r). (2.3) 



q=q 



We will discuss Ruijsenaars-type models associated with the classical root systems, 
namely the A^-i, Bj., Cr, D^. and BCr- The fact that all the roots of the classical root 
systems are neatly expressed in terms of the orthonormal basis of W makes formulation 
much simpler than those systems based on the exceptional root systems. The sets of posi- 
tive roots of the classical root systems are 



Ar-l 


A+ = 


= {cj - Cfcll < i < < r}, 








Al+ 


= {cj ±ek\l<]<k< r}, 


A5+ = 


{ej\l <]< r}, 


Cf 


A5+ 


= {Cj ±ek\l<i<k< r}, 


A^+ = 


{2e,\l< j<r} 


Dr 


A+ = 


= {cj ± efc|l < j <k<r], 








AAf+ 


= {Cj ± Cfell <j<k<r} 








A5+ 


= {e,|l <J <r}, Al+ = 


{2e,|l < 


j < r}, 



where {ej} is an orthonormal basis of W. The subscripts L, M and S stand for long, middle 
and short roots, respectively. 

It is well known that the non-simply laced root systems are obtained from simply laced 
ones by folding. In the present case, the relevant ones are: 

^2r-l — *■ Cr, Dr-\-l — ^ Bj-, — ^ BCr- (2-4) 

At the level of the dynamical variables and Hamiltonians, the above foldings are realised as: 

(j = l,...,r), (2.5) 

(2.6) 

-^j, (j = l,---,?^), Pr+l = gr+l = 0.(2.7) 

2.1 Calogero and Sutherland Systems 

For later comparison, we summarise the Calogero and Sutherland systems associated with 
the classical root systems only, namely the i?,., C,., and BCr- 

The Hamiltonian of the classical Calogero and Sutherland systems is 



^2r-l 




P2T + 1- 


-j ~ Pj^ 


Q2r+l~j 


Dr+l 


^Br 


Pr+1 -- 


= Qr+l = 0, 




A2r ^ 


BCr 


P2r+2- 


-J ~ ~Pj^ 


<l2r+2-j 



1 

^cs(p,g) = ^$^P? + K^(g), (2.^ 



2 

i=i 

where the classical potential Vq is given below explicitly. For all cases this classical potential 
Vc can be written in terms of the prepotential W{q) ^7] 

The equations for the equilibrium positions ()2.3|) reduce to the following equations: 

dW{q) 



2.1.1 Calogero Systems 



= 0, (j = l,...,r). (2.10) 

q=q 



The classical potential Vq and prepotential W are 

,2 1 «2„2 



j=i peA+ 



r 



^(^) = -2^^!+ (2.12) 

j=i peA+ 



where uo is the (positive) frequency of the harmonic confining potential and Qp are real 
positive coupling constants depending on the length of the roots. They are: one coupling 
g for all roots for the A^-i and D^.^ two independent couplings and gs for Bj. and Cr 
corresponding to the long and short roots, respectively. These conventions are the same for 
all other types of potentials considered in this paper. For example, the Cr model is 

= -| E ^1 + f E log k| - +9Lf: log |2g,| . 

j=i j,kji j=i 

There is no distinction between the rational B^. and Cr models because of the replacement 
gs ^ Ql- The Dr model can be considered as a special case of the Br with gi = 9 and 

The systems obtained by folding (j2.5|) - ()2.7p have a special ratio of couplings. They are 



folded : {gL.gs) = {\,l)g, folded 5^ : {gL.gs) = {l,2)g, 

folded BCr : {gi, 9m, 9s) = 1, l)fl'- 

The equations for the equilibrium position ()2.10|) are 



(2.13) 



fe=i 



Ar-i : y^-_ - = -gi, (2.14) 

Br : (2.15) 

"t^Qj-Ql 9l ' 9l qj 

Cr : E^ = -g.-^-, (2.16) 

iriq]-ql 9s 9s 

k^j 

Dr : j:-^,=--q,- (2.17) 



They determine the zeros of the Hermite and Laguerre polynomials. In other words, if 
we define = \f^yj for A^-i, then the polynomial having {yj\ as zeros is the Hermite 
polynomial [1011111111: 

2' ri(x - ») = r! t ^^^^^ S H„(,). (2,18) 



For the Cr (or Br) model let us define qj = \/^yj, ^ ~ ^ ~ ^ then having {yj} as zeros is 
the Laguerre polynomial (TUl dl E] : 

fri(-«i)-i:(::;)^s4<..(,). (2.19) 



j=l j=0 

For the D^. root system, it is the Laguerre polynomial Li~^\x). 
The identities between the Hermite and Laguerre polynomials 

2-^'H2r{x) = (-l)V!L!r^\x2), (2.20) 
2-^'-'H2r+i{x) = xi-lYr\Li^\x^), (2.21) 

are well-known. The former identity ()2.20|) for the even degree Hermite polynomial can 
be understood as a consequence of the folding of the root system ^2^-1 — > Cr, see fl2.5|) . 
Likewise the latter identity (j2.2H) for the odd degree Hermite polynomial can be understood 
as a consequence of the folding of the root system BCr, see (|2.7|) . Next let us consider 

the folding -Df+i Br (j2.6j) . which leads to the identity [4 among the Laguerre polynomials 
of different indices: 

{r + l)Li-'^{x) = -xL^;\x). (2.22) 

We will see that these identities (j2.2m) . (j2.2H) and (j2.22p are also nicely deformed with one 
parameter (jOU)) . (jOg|) and (jOUj) and with two parameters (|07j) . and P3T|) . 

2.1.2 Sutherland Systems 

The classical potential Vq and prepotential W are (except for Vq of BCr) 

Vciq) = I J2 -4^, ^ (2-23) 
^^(g) = I] ^plog|sin(p-g)|, (2.24) 

PGA+ 

where (?p are real positive couphng constants. The classical potential Vc of the BCr model 
is given by 



2 f^, Vsin {qj - qu) sin (g^ + qu) 



r 



^ sm 20, 2 sm o,- 



The Br (Cr) potential is obtained by setting = 0, qm — ^ dL, ids = 0, Qm ds) 
The equations for the equihbrium position ()2.10|) are 



Ar-l 


fc=i 

r 


- Qk) 


= 0, 


Br 


■ $^(cot(g, 

fc— 1 
r 


-Qk) 


+ cot(gj + 


Cy" '■ 


: ^(^cot(gj 

fc=i 

kT^J 

r 


- qk) 


+ cot(gj + qk)j 




■ $^(cotfe 

fc=i 

kj^j 

r 


- Qk) 


+ cot(gj + qk)^ 


Dr 


fc=i 


- Qk) 


+ cot(gj + qk)^ 



(2.26) 



-— cotgj, (2.27) 
9l 



9s 



9m 9m 



(2.28) 



(2.30) 



The equihbrium positions of the Ar-i model are related to the Chebyshev polynomial and 
those of the other models are related to the Jacobi polynomials. 

For the A^-i, the equilibrium positions are equally- spaced" and translational invariant, 

g= - *(r,r- 1,---,1) + ^ *(1,1,---,1), ^ G M : arbitrary. (2.31) 

r 

r 

We choose this constant shift ^ such that the "center of mass" coordinate vanishes, qj = 0: 

* - ^^i^ - ^ - f - - 0- 1. .... (2.32) 

Then the degree r (the dimension of the vector representation) polynomial in x, having zeros 
at {singj}, 

T"' f[{x - sin qj) = fl(x- cos j 4£,f t,(x) , (2.33) 

j=i j=i ^ 

is the Chebyshev polynomial of the first kind, Tn{cosip) = coslrup). 

For the solution {qj} of the BCr fj2.29p . cos2gj is the zero of the Jacobi polynomial 
p!^"''^\x) with a = ^ + ^- land/5=^-l, 

3=1 j=0 
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It is easily shown that = f ~ Qj satisfies ()2.29j) with a ^ {3, which imphes Pn' (—x) = 
{—!)"' P!f'°'\x). For the solution {qj} of the Cr ()2.28|1 . cos2gj is the zero of the Gegenbauer 
polynomial cl°'^'^\x) with a = ^ — I, 

j=l j=0 

(2.35) 

This is a special case of Pr"''^\x) with another normalisation, 

C^"+^\x) = (2"+'^) ("+'^)"'pj"'")(a;). (2.36) 

For the solution {qj} of the 5^ ()2.27p . cos2qj is the zero of Pr-°'"^Ha;) with a = 2^ - 1. For 
the solution {qj} of the Dr ()2.30|) . cos2qj is the zero of Pr ^' ^''(a;). 

The known identities between the Chebyshev and Jacobi polynomials and between the 
Jacobi polynomials can be understood as consequences of the folding: 

A,r.,^Cr : 2'~'^T,r{x) = {-ir{'^;yPy'^'-'-\l-2x'), (2.37) 
A2r^BCr : 2-'^T2r+i{x) = x{-lY{l'y'py^'~'^\l-2x'), (2.38) 
A+1 ^ 5, : 2(r + l)Pi-l'-'\x) = r{x - l)Py-'\x). (2.39) 

We will see in the following that the first two identities are not deformed but the third one 
is nicely deformed ()5.53|) . 

The Gegenbauer polynomial and the Jacobi polynomial are also related by the quadratic 
transformations : 

Ct'\x) = 22"rH")(^:)-Vi"-^\2a;^-l), (2.40) 

ctj\x) = 2-+^r4r)(":T'^^i""^(2^' - 1) • (2.41) 

However these identities do not seem to have a folding type explanation. Indeed the deformed 
Gegenbauer polynomials ci^^^^x, 5) ()5.18p . C^~^'^\x,6) ()5.26|) and the deformed Jacobi 
polynomial Pn°''^\x,6) ()5.2|) do not satisfy this type of identities for generic 6. 

2.2 Ruijsenaars-type Systems 

Following Ruijsenaars-Schneider (Hj and van Diejen (?], the Hamiltonian of the Ruijsenaars 
systems is (the deformation parameter /3' of e^^'^ is set to unity, /?' = !) 

H{p,q) = i2 Uoshp,^V,{q)V;{q)-^{v,{q) + V*{q))) . (2.42) 
i=i ^ ^ 

9 



The form of the function Vj = Vj{q) and its complex conjugate V* are determined by the 
root system A as : 

r 

Ar-i : Vj{q)=w{qj)l[v{qj~qk), (j = l,...,r), (2.43) 

r 

Br, Cr, Dr, BCr : Vj{q) = w{qj) Y[v{qj - qu) v{qj + g^) , (j = 1, . . . , r). (2.44) 

fc=i 

The elementary potential functions v and w depend on the nature of interactions (rational, 
trigonometric, etc) and the root system A. Their explicit forms will be given below. When 
V satisfies the simple identity J2j {^jiQ) + ^/{q)) = const. , the Hamiltonian ()2.42p could 
be replaced by a simpler one 

r 

H'{p,q) = J^coshp, ^V,{q)V*iq), (2.45) 
i=i 

which is obviously positive definite and usually used as a starting point for the trigonometric 
(hyperbolic) interaction theory. 

The above Hamiltonian ()2.42j) is a hyperbolic function of the momentum operator p 
instead of a polynomial in the hierarchy of C-M systems or other ordinary dynamical systems. 
In quantum theoretical setting this Hamiltonian causes finite shifts of the wavefunction in 
the imaginary direction, ie coshp?/'(g)= ^{iplq — ih) + iplq + ih)), in which h is Planck's 
constant. This is why the R-S systems are sometimes called 'discrete' dynamical systems. 

The equation ()2.3|) of equilibrium positions ()2.2p can be simplified in the following way. 
By expanding the Hamiltonian around the stationary solution ()2.2|) . we obtain 

1 

H{p, q) = K{p) + P{q) + higher order terms in j9, K{p) = - N I Vj (g) | , (2.46) 

i=i 

and the 'potential' P is given by 

P{q) = E (^vMv;{q) - \ [vM) + v;{q))) =-\Y. (\A^ " \/^)' • (2-47) 

This should be compared with the classical potential in the Calogero- Sutherland systems 
(12. 9p . It is obvious that the equilibrium is achieved at the point (s) in which all the functions 
Vj become real and positive: 

= ^*(?)>0, (j = l,2,...,r). (2.48) 
10 



The equilibrium point is the absolute minimum of the potential P. The system of equations 
fl2.48|) is invariant under any permutation of {(jj}- For v and w considered in this paper, 
they are also invariant under the transformation q ^ q' = —q. Except for the A case, they 
are also invariant under the transformation q q' = ^{^iqi, ■ ■ ■ , ^rqr) with = ±1. 
The functions v and w considered in this paper have properties 

v{—x) = v*{x), w{—x) = w*{x), (2.49) 

which allow the folding (j2.5p - (|2.7|) of the Ruijsenaars type Hamiltonians: 

w{x) =w^{x)v^{2x), (2.50) 
w{x) =w^{x)v^{x)v'^{2x),{2.51) 
w(x) = w°{x)v^{xY, (2.52) 

where H is 

H{p, q) = G^shp, ^Jv,{q)V*{q) - ]- {v,{q) + V;{q)) 
j=i V 

r 

Vj{q) = w{qj) Yl v{qj - qu) v{qj + qt) . (2.53) 

fc=i 

The folded systems fl2.50p . ()2.5H) and ()2.52|) correspond to the folding A2r-i — > Cr, — ^ 
BCr and -Dr+i Bj., respectively. The coupling constants in these folded systems are not 
independent as shown in ()2.13p . 

2.2.1 Ruijsenaars-Calogero Systems 

The first example to be discussed is the 'discrete' analogue of the Calogero systems PP, 
to be called the Ruijsenaars-Calogero systems, which were introduced by van Diejen for 
the classical root systems only [71 |H] . The original Calogero systems pOj have the rational 
(l/(distance)^) potential plus the harmonic confining potential, having two coupling con- 
stants g and uj for the systems based on the simply-laced root systems, A and D, and three 
couplings u and gi for the long roots and gs for the short roots in the B and C root systems. 

Two varieties (deformation) of 'discrete' Calogero systems are known. The first has two 
(three for the non-simply-laced root systems) coupling constants g {gi and gs) and a which 



H^^^-'{p,q) 
H^'^{p,q) 
H''^^\p,q) 



12r+l- 



P2r+2-j=-Pj 
12r+2-j=-1j 



Pr+l=0,qr+l=0 



2H{p,q), 



2H{p,q), 



Hip,q), 



ViX) = V (X) 



V[X) = V [X), 



v(x) = v^(x) 
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corresponds to uj in the Calogero systems. The second has three (four for the non-simply- 
laced root systems) coupling constants g [ql and gs) and a, b both of which correspond to 
u. The integrability (classical and quantum) of these systems was discussed by van Diejen 
in some detail [HIH]. The former can be considered as a limiting case {b cxo) of the latter. 

Linear Confining Potential Case 

The dynamical system is defined by giving the explicit forms of the elementary potential 
functions v and w. For the simply-laced root systems A and D they are: 

A, D : v(x) = l-i — , w(x) = l + i — , (2.54) 

X a 

in which a and g are real positive coupling constants. For the non-simply-laced root systems 
B, C and BC, we have: 

w{x) = (l + i-)(l-i^£), (2.55) 



B : 


: v[ 


[x) 


= 1 


-i — 

X 


C : 


: v[ 


[x) 


= 1 


-i — 

X 


BC : 


: v[ 


[x) 


= 1 


-i — 








X 



a 



w(x) = + , (2.56) 

^x)=(l + .|)(l-.|)(l-z|), (2.57) 

in which a, gi, gs, do, Oi, 92 are independent real positive coupling constants. Normalisation 
of the coupling constants is chosen such that they reduce to those of the Calogero models in 
the small coupling limits discussed below. The D model is obtained from the B model by 
gL = g and gs = 0. The B and C models are special cases of the BC model. In contrast 
to the Calogero case, those based on the B and C systems are different. The difference of 
the length of the roots is immaterial, since it can be absorbed by the coupling constants 
normalisation. It is rather elementary to derive the forms of the Hamiltonians of the B and 
C systems from those of the D and A systems by folding (see (|2.52p and (|2.5(J|) ). In all these 
cases the 'potential' P ()2.47j) grows linearly in |g| as |g| — oo. Except for the BC case, there 
are simple identities : J2j i^ji^l) + — const. 

In the limit of small coupling constants, namely, by recovering the deformation parameter 

{Pj, -, 9, 9l, 9s, 9o, 9i, 92) {Pj, -, 9, gL, 9s, 9o, 9i, 92), (2.58) 
and taking /5' — > limit, the Hamiltonian ()2.42j) tends to that of the corresponding classical 
Calogero system ()2.8|1 with = \ {BCr tends to Cr with gs = go, 9l = 9i+ 92) 

■^H{p, q) i/calogero(P, ?) + COUSt. (2.59) 
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It is interesting to note that the equations determining the equihbrium ()2.48|) . in general, 
can be cast in a form which looks similar to the Bethe ansatz equation. For the elementary 
potential ()2.54|) - ()2.57|) . the equilibrium positions {qj} are determined by: 

A,_, : fj ^-^^-W = (2.60) 
f-J^ Qj -Qk + W a + iQj 

Qj -Qk- igL qj + qk-igL a- iqj f 2qj + igs 



. TT qj - qk- i-yL qj -r qk- i-ul ^ » - i-qj i ^qj ^ tys \ ^2 g]^^ 

IJ- Qj -qk + igi qj + qk + igL a + iqj \2qj - igs) ' 



^ . TT Qj -qk- igs Qj + qk- igs ^ « - iqj qj + igL 
IJ- qj -qk + igs qj + qk + igs a + iqj qj -igi' 



. TT qj - qk- igo qj + qk- igo ^ a - iqj qj + igi q-j + 
.1 Qj -qk + igo qj + qk + igo a + iqj qj - igi qj -ig2' 



. -TT qj -qk-ig Qj + qk-ig ^ a - iqj 

■JJ- qj -qk + ig qj + qk + ig a + iqj ' 

In the above small coupling limits, these equations reduce to ()2.14|) - ()2.17p . Thus the Bethe 
ansatz-like equations ()2.60|) - ()2.64|) would give deformation of the Hermite and Laguerre 
polynomials, as we will see in section ITT] and I01 They are not the so-called q- deformed 
Hermite or Laguerre polynomials 



Quadratic Confining Potential Case 

In this case the elementary potential function v is the same as before, but w is different. For 
the simply-laced root systems A and D, the elementary potential functions are: 

A, D : v{x) = l-i^, w{x) = + + {a, b, g > 0). (2.65) 

For the non-simply-laced root systems B, C and BC, we have {gi, gs, go, gi, g2 > 0): 



B : 


: v{x) 


= 1 


-i — 

X 


C ; 


: v{x) 


= 1 


-i — 






X 


BC ; 


: v{x) 


= 1 


-i — 






X ' 



a 



l + i — + , (2.66) 



2 



b J\ 2xJ 



'-7)(^-T)(-'f)(-'f)- 

The D model can be considered as a special case of the B model hj gi = g and gs = 0. 
The B and C models are special cases of BC model. As in the previous case, the forms of 
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the elementary potential function w for B and C systems are determined from those of the 
D and A systems by folding ()2.52|) . ()2.50p . In all these cases the 'potential' P ()2.47|) grows 
quadratically in |g| as |g| — oo. The small coupling limit ()2.58|) (and | — ^ y) gives the same 
classical Calogero systems as before ()2.59|) . ()2.8|) with uj = ^ + ^• 

The equations ()2.48|) determining the equilibrium positions {qj} for the elementary po- 
tential ()2.65|) - ()2.68|) are expressed in a form similar to the Bethe ansatz equation: 

_"7 Qj - Qk + ig a + iqj b + iqj 



^ . TT Qj -Qk- iQl Qj + qk- Wl ^ a - iqj b - iqj f 2qj + igs \ 

M Qj -qk + igi qj + qk + igL a + iqj b + iqj \2qj - igs J ' 



^ . TT Qj -Qk- igs Qj + qk- igs ^ a - iQj b - iqj qj + igi 
Qj -qk + igs Qj + qk + igs a + iqj b + iqj qj -igi' 



r 



. TT Ij -Qk- igo qj + qk- igo ^ a - iqj b - iqj qj + igi q^ + ig2 ,^ 
fej- Qj -qk + igo qj + % + igo « + iqj b + iqj qj - igi qj -ig2' 

^ . TT Qj -Qk-ig qj + qk-ig ^ a - iqj b - iqj 
■Jli qj -qk + ig qj + qk + ig a + iqj b + iqj ' 

k^j 

They define another type of deformation of the Hermite and Laguerre polynomials, since the 
small coupling limit of the above Bethe ansatz-like equations gives the same equations as 
before (j2.14p - (j2.17p . determining the zeros of the Hermite and Laguerre polynomials, with 
u = ^ + I- These will be discussed in sections U.2\ and 

2.2.2 Ruijsenaars-Sutherland Systems 

The discrete analogue of the Sutherland systems 0, to be called the Ruijsenaars-Sutherland 
systems, was introduced originally by Ruijsenaars and Schneider jHj for the A type root 
system. The quantum eigenf unctions of the A type Ruijsenaars-Sutherland systems are 
called Macdonald polynomials ^H] , which are a one-parameter deformation (g-deformation) 
of the Jack polynomials Here we will discuss the Ruijsenaars-Sutherland systems for all 
the classical root systems. A, B, C, D and BC 8J. The structure of the functions {Vj{q)}, 
()2.4H|1 and ()2.44j) are the same as in the Ruijsenaars-Calogero systems, but the elementary 
potential functions v and w are trigonometric instead of rational. Because of the identity 
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Y7j=i i^jil) + ~ const., the Hamiltonian ()2.42j) could be replaced by a simpler one 

The elementary potential functions v and w are: 



A, D : 


: v{ 


» 


= 1 


— ztanhgfcota;, 


w{x) ■■ 


= 1, 




(2.74) 


B : 


: v{ 


[x) 


= 1 


— i tanh ql cot x, 


w{x) 


= (1 


— i tanh ^ cot x)^, 


(2.75) 


B' : 


: f 1 


[x) 


= 1 


— i tanhgf^ cot x, 


w{x) 


= 1 - 


- i tanh cot x, 


(2.76) 


C : 


: t;( 


[x) 


= 1 


— i tanh5f5 cot x, 


w{x) 


= 1 - 


- ita.nh.2gL cot 2x, 


(2.77) 


c 


: v{ 


[x) 


= 1 


— i tanh5f5 cot x, 


w{x) 


= (1 


— i tanh (7/, cot 2x)^, 


(2.78) 


B'C ; 


: t;( 


[x) 


= 1 


— i tanh (7m cot x, 











w{x) = {I — i tanh (75 cot x)(l — i tanh 2(7 j;^ cot 2x), (2.79) 

with similar coupling constant notation as in the rational cases. Normalisation of the coupling 
constants is chosen such that they reduce to those of the Sutherland models in the small 
coupling limits discussed below. The D model is the special case of the B model hj ql = g 
and Qs = 0. The B' and C models are special cases of the B'C model. As in the rational cases, 
the forms of the elementary potential function w for the B and C systems are determined 
from those of the D and A systems by folding (j2.52j) . (j2.5U|l . 

The original Sutherland models are obtained in the limit in which all the coupling con- 
stant (s) become infinitesimally small. By recovering the deformation parameter (3\ 

{Pj, 9, 9l, 9m, 9s) {Pj, 9: 9l: 9m, 9s), (2.80) 

and taking — > limit, the Hamiltonian ()2.42|) tends to that of the corresponding classical 
Sutherland system ()2.8|) 

■^H{p, q) -f/'suthorland(P, Q) + COUSt. (2.81) 

In 'strong' coupling limits, g, gi, gM, 9s +00, the elementary potential functions v 
and w take simple forms: 

v{x) ^ 1 — icotx, w{x) ^ 1, 1 — icotx, (1 — icotx)^, etc. (2.82) 

The deformed polynomials take simple forms in this limit as we will see in section El 
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The equations ()2.48|) determining the equihbrium positions {qj} for the elementary po- 
tential ()2.74p - ()2.79p are expressed in a form similar to the Bethe ansatz equation: 

TT tan((7,- — Qh) — i tanh q 

Ar-i : n , L> / = ^' 2.83 

tan(gj — g^j + ^tanh^f 



^ . -TT tan(gj - qk) - itanh^^ tan(gj + g^) -itanhgL _ f tan q^ + i tanh — x 
^ "{li tan(gj — qk) + i tanh tan(gj + g^) + ^ tanh (7^ \tan — i tanh ^ ' ' 



, ^ tan{q j - qk) -i tanh ql tan{q j + qk) - i tanh ql tan g^ + i tanh 5(5 /o ocr\ 

-D^ : II 7^ — : T-z : : = z : \ , [Z.ooj 

tan(gj — qk) + i tanh tan(gj + qk) + i tanh (yf^ tan qj — i tanh gs 



kltj 



^ . -TT tan(gj - g^) - itanhgfg tan(gj + qk) - itanh^g _ tan 2gj + i tanh 2gL ,^ 
^ "Jli tan(gj — (jk) + i tanh gs tan(gj + g^) + i tanh gs tan 2gj — i tanh 2gL ' 



^, . -TT tan(gj - g^) - itanhgfg tan(gj + g^) -itanh^g ^ / tan2gj + itanhgfz, V^^ 
^ _ tan(gj — (jk) + i tanh (jfg tan(gj + (jk) + i tanh (yfg \tan 2gj — i tanh (7^ / ' 



kT^i 



^'^ . tan(gj - qk) - i tanh 5(4/ tan(gj + g^) - i tanh 



^ tan(gj — g^) + itanh5fj\/ tan(gj + qk) + itanhfyf^ 

tan (jj + i tanh gs tan 2gj + i tanh 2gL 
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tangj — itanh^fg tan2gj — ita,\^\i2gi 

JJ . -TT tan(gj - g^) - itanhgf tan(gj + qk) - itanh^ ^ ^ 
^ -L-L tan(gj — g^) + itanh^f tan(gj + g^) + i tanh (7 



From the property mentioned after (j2.48j) and the fact that (j2.83p - (j2.89p are the equations 
of {tangj}, we can restrict qj to < qj < it/ 2 (except for the A^-i case). In the small 
coupling limit, these equations (j2.83j) - (j2.89j) tend to (j2.26j) - (j2.3(Jj) . Thus the Bethe ansatz- 
like equations (|2.83|) - (|2.89|) would give deformation of the Chebyshev and Jacobi polynomials 
as we will see in section El 

2.3 Orthogonal Polynomials and Three Term Recurrences 

It is well known that orthogonal polynomials satisfy three-term recurrence ^Hl 120] and con- 
versely a sequence of polynomials satisfying three-term recurrence are orthogonal with re- 
spect to certain inner product with some weight function. Here we will introduce appropriate 
notation by taking the classical orthogonal polynomials as examples. 

Let {/n(a^)}J^Lo be a sequence of orthogonal polynomials with real coefficients. That is 
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is a degree n polynomial in x and they are mutually orthogonal {fn, fm) = hnSn,m 
{hn > 0) with respect to an (positive definite) inner product {f,g) = J f{x)g{x)w{x)dx 
(w(x) is a weight function). Let /"""'^(x) be a monic one, /„(x) = c„/™°°"^(x) = c„(x" + - ■ ■). 
Then /™°°"^(x) satisfies three-term recurrence : 

/~(x) - (x - a„)/— (x) + &„/rr (x) = , (n > 0) , (2.90) 

where we have set /™°°"^(x) = ^, and a„ (n > 0) and 6„ (?t- > 1, &o is unnecessary, 6„ > 0) 
are real numbers. The constants a„, 6„ and hn are given by 

(a;/„(x),/„(x)) , cl_^ (/„(x),/„(x)) 



/i„ = (l,l)4j]6,. (2.91) 



(/n(x),/„(x)) ' C2 (/„_l(x),/„_i(x)) .^^ 



Namely fn{x) satisfies the three-term recurrence 



— /„+i(x) -(x-a„)/„(x) + 6„^/„_i(x) = 0, (n>0). (2.92) 



Cra+l C^— 1 

For a„ = (n > 0) case, /n(x) has a definite parity, /n(— x) = (— l)"'/„(x), and the constant 
term of the even polynomial is 

n 

/2n(0) = (-l)"C2„n^2,-l. (2.93) 

i=i 

Conversely, if {/n(x)} is defined by the three-term recurrence (j2.9U|) . namely, real numbers 
a„ (n > 0), bn {n > 1, > 0) and c„ (n > 0, c„ 7^ 0) are given, then {/n(x)} is a sequence 
of orthogonal polynomials with respect to some (positive definite) inner product (■,■)• 

In the rest of this subsection we summarise the three-term recurrence, generating func- 
tions, differential equations, etc for the Hermite, Laguerre and Jacobi polynomials for later 
comparison with the corresponding quantities of the deformed polynomials. 

The Hermite polynomials Hn{x) (j2.18|) are orthogonal with respect to the inner product 
(/)fl') = /(^)5'(2^)6~^^'^3;, {hn = 2"'n!y7r) and satisfy the three-term recurrence ()2.9()j) 
with 

Ti 

a„ = 0, bn = -, c„ = 2", (2.94) 
Hn+iix) - 2xHn{x) + 2nHn-i{x) = 0. (2.95) 



^Hereafter we adopt the convention ~ and /o(a;) = 1 for all the polynomials in this paper. 
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The generating function and orthogonality are 

G{t, x) = J2 -Hn{^) = e-*'+'"*, {G{t, x), G{s, x)) = (2.96) 



n=0 



The Laguerre polynomials llC^ [x) ()2.19|1 are orthogonal with respect to the inner product 
{.f-,9) = f{x)g{x)x"'e~^dx, {hn{a) = T{a + n + l)/n\ , Rea > —1) and satisfy the three- 
term recurrence ()2.9()j) with 

an = 2n + l + a, 6„ = n(n + a), c„ = (-l)"/n!, (2.97) 
(n + l)4+i(a;) + {x-{2n + a + l))Li^\x) + (n + a)L^^\{x) = 0. (2.98) 

The Jacobi polynomials Pn"''^\x) ()2.34p are orthogonal with respect to (/,(?) = J^^ f{x) 
g{x){l - x)"(l + x)'^dx, {hnia,(3) = ^^^^^I+i ^'^"r^rit+p+TXiT^ ) satisfy the three-term 
recurrence ()2.90|) with 



, 6„ = 4n(n + «)(n + /?) X , (2.99) 

«2n«2ri+2 "2ra-l "2?! "2 



^2n "2n+l 

c„ = 2-"("+'3„+2"), d^ = a + P + m, (2.100) 
2(n + 1) (« + /? + n + 1)(« + /5 + 2n)Pi°f (x) 

-(a + /3 + 2r2 + 1) ({a + P + 2n){a + (3 + 2n + 2)x + - /J^) P^^'^^ (x) 

+2{a + n){l3 + n){a + (3 + 2n + 2)Pi°f (x) = . (2.101) 

The Gegenbauer polynomial ci^^^^x) ()2.35|) is a special case of the Jacobi polynomial 
C^"+^\x) = (^"+")("+")"^Pi"'"Hx). The Chebyshev polynomial of the first kind T„(x) is 
also a special case of the Jacobi polynomial T„(x) = 2~^(^"^^) ^pi ^' ^''(x). The Legendre 
polynomial is P„(x) = Pn'^\x) 

The differential equations of the Hermite, Laguerre and Jacobi polynomials are 

^Hn{x) - 2x±Hn{x) + 2nHn{x) = , (2.102) 

x^4") (x) + (a + 1 - x) (x) + (x) = , (2.103) 
(1 - ^')&Pi''''\^) - (« + /? + 2)x) f Pi"-'^)(x) 

+n(n + a + /? + l)Pi"''^)(a;) =0. (2.104) 
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3 Deformation of the Hermite Polynomial 

3.1 Linear Confining Potential Case (1 Parameter Deformation) 

For the solution {qj} of the Ar-i equation ()2.60|) . let us define 

and introduce a degree r polynomial in x having zeros at {yj}: 

r 

Hr{x,S)''^rl[{x-yj). (3.2) 

It is a deformation of the Hermite polynomial (|2.18p such that 

lim HJx, 5) = HJx). (3.3) 

(5^0 

If {qj} is a solution of ()2.60|) . so is {—qj}, which would imply that the deformed polynomial 
Hr{x, 6) has a definite parity 

Hr{-x,6) = {-iyHr{x,6), (3.4) 

as with the original Hermite polynomial Hr{~x) = {—lYHr{x). 

The equation for the equilibrium ()2.60|) can be written as (we replace r by n) 



n 



Tjyi^yh^ = y_ill^, (3.5) 

From this equation, we obtain the following functional equation for Hn{x, 6) {n > 1), 

(^x + i-j^Hn{x + iv^, 5) — — i—j^Hn{x — iVS, 6) = 2iAnHn{x, 6), (3.6) 

because the LHS is i times a degree n polynomial in x with real coefficients which vanishes 
at X = yj. Here An = An{S) is a real constant. This functional equation contains all the 
information of the equilibrium. The number of unknown coefficients (coefficient of x^ term 
of Hn{x, 6) {k = 0,1, . . . ,n — 1) and An) and the number of equations (coefficient of x^ term 
of ()3.6|) {k = 0,1, . . . ,n)) are both n. The constant An is determined by the coefficient of x"^ 
term of this equation, 

An = ^{l + nS). (3.7) 
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The functional equation ()3.6|) can be written as a difference equation, 

^',i.v^^-(^''^) - '^^D.,^V5Hn{x,5) + 2nH^{x,5) = 0, (3.8) 
where D^^h is a (central) difference operator, 

g../W = ^'" + ";/'""" . (3.9) 
In the 5^0 limit, (j3.6p reduces to the differential equation of the Hermite polynomial 

(inna). 

The uniqueness (up to normalisation) of the solution of the functional equation (j3.6j) 
is easily shown (Proposition IA.1|) . Therefore it is sufficient to construct one solution of 
(j3.6p explicitly. This is done by using the three-term recurrence (see Appendix). The result 
is as follows; The functional equation (j3.6|) implies that the deformed Hermite polynomial 
Hn{x,6) satisfies the three-term recurrence (j2.9(J|) with 

71 / Tl — 1 \ 

a„ = 0, 6„ = - 1 + ^—5 , c„ = 2", (3.10) 



2 V 2 

Hn+i{x,5)-2xHn{x,5) + {2n + n{n-l)5^Hn-i{x,5) = 0. (3.11) 

Since 5 is positive in this case ()3.1|) . 6„ is also positive. Therefore Hn{x,6) is a set of 
orthogonal polynomials with respect to some positive definite inner product. 

Here we present another derivation of this three-term recurrence. Let us consider the 
generating function 

°° t" 

G{t,x,6) = y^-H4x,S), (3.12) 

n=0 

which satisfies 

({l + 5t^)— + 2{t-x)^G{t,x,6) = 0, (3.13) 

as a consequence of the three-term recursion (j3.11|) . This linear differential equation with 
the initial condition G{0, x,6) = 1 can be easily solved and we obtain 



exp(2x 



arctan \/6 1 ^ 



G{t,x,5)= ^f. . (3.14) 

(1 + dt^)s 

In the (5 — > limit, this generating function tends to that of the Hermite polynomial ()2.96|) . 

lim G{t, X, 5) = e-*'e2"* = 6-*'+^-* = G{t, x). (3.15) 

<5— >0 
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The functional equation of G(t,x,6) is obtained from ()3.6p : 

(x + i^^G{t,x + iV6,6)- (^x-i^^G{t,x-iV6,6) = ^(^1 + 6t^^G{t,x,S). (3.16) 

Since the solution of (j3.6|) is unique (up to normalization), it is sufficient to show that (j3.14j) 
satisfies this functional equation. This can be easily verified by explicit calculation, in which 
the following formula derived from ()3.14|1 is useful, 

G{t, X ±1^/6,6) = ^^'^^ Git, X, 5). (3.17) 
1 =F iVSt 

The explicit form of the inner product {f,g) = f{x)g{x)w{x,6)dx , ie the weight 
function w{x,6) is determined by using the generating function. Here we list main results 
only without derivation. Let us fix its normalisation by (1, 1)5 = a/vt- From the general 
theory (j2.9H) . the orthogonality of Hn{x,6) is 

n-l 

{Hn{x, 6), Hm{x, 5))s = (5„,„/i„, K = ^Tn\ J](l + ij5), (3.18) 

which leads to 

(G'(t,a;,(5),G(s,x,(5))5 = (1 - (5ts)-i. (3.19) 
Here we have used the identity 

00 ^ n—\ 

^(^) = E ^ 11(1 + \3b) = (1 - \8x)--^. (3.20) 

n=0 ■ j=0 

The weight function is expressed as 

/ /"°° cos2xt 2i"^^/l X 1 x\ , ^ 

w(x,(5) = / ———-^ = —=b[-^ I— .--%—] (3.21) 



v^(coshv^t)^ ^5 5 a/5 



2t-^r(i + .^)r(i-.^) 2l-Mr(i + .:^ 



X \ \2 



r(|) r(|) 



(3.22) 



The undeformed limit of the weight function lim^^o w(a;, 5) = e ^ can be verified by using 
the asymptotic expansion of the F- function. The Taylor series of w(x, 5) in powers of 5 reads 

w(x, b) = e"^' ( 1 + — (3 - 12x2 + 4x^) + (45 - 1320^^ + 2280x^ - 864x^ + 80x^) 

^ ' ^ V 24^ ' 5760^ ^ 

+ (-14175 - 71820x2 + 865620x^ - 1042272x^ + 386928x^ 

2903040^ 

-52416x1° + 2240x^2) + . . . j . (3.23) 
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Among many interesting properties of Hn{x, S), we present only 



H2ni0,6) = {-l)^^l[{l+j5), see (Eli, (3.24) 



3=0 



(ii) ^H^ix, 5) = 2 ^ (2A:)! L ^ {-5f H^_,_,,{x, 5), (3.25) 
^ fc=0 ^ ^ 

which is a deformation of -^Hn{x) = 2nHn{x). 

Remark : We may take the three-term recurrence (|3.11|) as the definition of the deformed 
Hermite polynomial Hn{x,S) for an arbitrary (complex) parameter 6. Then Hn{x,6) is a 
polynomial in x of degree n and in 5 of degree [|] with integer coefficients. We will not 
repeat similar remarks which are valid for almost all of the deformed polynomials in this 
paper. 

3.2 Quadratic Confining Potential Case (2 Parameter Deforma- 
tion) 

For the solution {qj} of the Ar-i equation (j2.69j) . let us define 

Qj = Vagyj, S=-, £ = j-, (3.26) 
a 

and introduce a degree r polynomial in x having zeros at {yj}'- 

r 

Hr{x,S,e)''^2^l[{x-y,). (3.27) 

It is a further deformation of the deformed Hermite poljTiomial defined previously, 

\imHr{x,6,e) = Hr{x,6), Hr{-x,6,e) = (-l)"iJ,(x, 5, e). (3.28) 
The symmetry between the two parameters a 6 is expressed as 

Hr{x, Se, e-^) = e^^Hrie'^x, S, e). (3.29) 



If we define Hr{x, Si'= ^ = S, ^2 *= f = f ) == Vl + e'' Hr{ ^^_^_^ , S, e), then this symmetry i 
more manifest, Hr{x, 61,62) = Hr{x,62,6i). 

The equation for the equilibrium ()2.69|) can be written as (we replace r by n) 

fjyi^yt^ = J_il^y_il^_ (3.30) 
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From this equation, we obtain the following functional equation for Hn{x,6,e) (n > 1), 

= 2{AnX^ + BnX + Cn)Hn{x, 6, e) , (3.31) 

because the LHS is a degree n + 2 polynomial in x with real coefficients which vanishes at 
X = Uj. Here = ^4^(5, e), i?^ = e) and Cn = Cn(5, e) are real constants: 

v4„ = l, S„ = 0, C„ = -((^+n)£-^+n+|n(n-l)5) . (3.32) 

This functional equation contains all the information of the equilibrium. The above func- 
tional equation ()3.3H) reduces to that of Hn{x,6) ()3.6|) in a proper limit e —>■ 0. This 
functional equation can be written as a difference equation, 

il-6ex^)Dl,^^^Hnix,6,e)-2{l+E)xD^^^^,H^{x,6,E) + 2n{l + {l + ^6)E)H^ix,6,E) = 0. 

(3.33) 

The functional equation 1)3.311) implies (see Appendix) the three-term recurrence ()2.9()j) 
for the deformed Hermite polynomial Hn{x,6,E) with 

n / n — 1 \ / n — 1\ d„ 
an = 0, bn = -{l + ^—6)(l + ^—SE)- ^ , c„ = 2^ (3.34) 



2 V 2 / V 2 / d2n-ld2n+l 

d^ = l+[l + ^-6jE, (3.35) 

/ \ (1 + ^5E)dn 

Hn+iix,S,E) -2xHn{x,6,E) + 2n + n(n - 1)5 ' g„_i(x,5,£) = 0.(3.36) 

^ ^ «2n-l«2n+l 

Since 6 and e are positive in this case ()3.2(jp . 6„ is also positive. From this three-term 
recurrence we obtain the differential equation for the generating function G{t, x, S, e) = 

E—Hn{x,6,E), 
71,1 



n=0 



( 8 4: \ 

[{d2n-ld2n+l) (— ^ 2x) + {- d2n~ld2n+lbn) t ) G (t , X , 6, e) = , (3.37) 

V 8t n ^ 

which is a third order linear differential equation with respect to t. The special case 5 = 
gives the original Hermite polynomial, Hn{x, 0,e) = (vT+T )^"if„(vT+Tx). The value at 
the origin of the even polynomial shows a characteristic deformation pattern, see ()2.93j) : 

^2n 0,5,£ = -1 j- M ^ , , TT^. 3.38 

^- J=o 1 + (1 + (j + - 
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4 Deformation of the Laguerre Polynomial 

4.1 Linear Confining Potential Case (1 or 2 Parameter Deforma- 
tion) 

BCr • For the solution {qj} of the BCr equation ()2.63|) . let us define 

/ X 3o 91+92 ^ 9i92 

Qj = Vm)yj, d = —, a = 1, 7 = , (4.1) 

« 9o 9o 

and introduce a degree r polynomial in x, having zeros at 

^l"na^,7,5)^=^^n(^-y|). (4.2) 

It is a two parameter deformation of the Laguerre polynomial such that 

\imLi''\x,j,S) = Li^'^x). (4.3) 
The equation for the equilibrium ()2.63|) can be written as (we replace r by n) 

M (%• + - yl + - ^71 1/| - ^(« + l)^5y^-^5' 

From this equation, we obtain the following functional equation for Ln"''(x,7,5) {n > 1), 



y + 2^=) + ^(a + l)v^y - 75)4-) ((y + tVs f, 7, 5) 



.v^W .1 

2t{A^y' + B^)Ll^\y',^,S), (4.5) 



1/ + - (y' -z{a + l)V6y- 7^)4"^ ((y - ^V~6)^ 7, ^ 



because the LHS is i times a degree 2n + 2 even polynomial in y with real coefficients which 
vanishes at y = ±yj. Here An = Al^\'y, 6) and i?„ = B^\'y, 6) are real constants: 

I I v/n -i- n/ -i- - \n \ « = 

This functional equation contains all the information of the equilibrium. The functional 
equation can be written as a difference equation, 



A„ = ^ (1 + (2n + a + i)5) , 5„ = ^ (a + 1 - 27 + (n + 7)5) • (4.6) 



((! + (« + + - 7)5 + |75^)^J,l,^LW(y^7, 5) 

+ (-2y3 + (2a + 1 + (27 - a - 1)5 )y + 75 ^ ^..^/^ 4"^ 7, S) 

+ iAny' + n5)L^:\y',^,5) = 0. (4.7) 
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In the (5 —i> limit, it becomes, 

y^,L^:\y') + (2a + 1 - 2y')^Ll:\y') + AnyL^:\y^) = 0, (4.8) 

which is equivalent to the differential equation of the Laguerre polynomial ()2.1()H|1 . 

By the same argument as Hn{x,6), see Appendix, the functional equation ()4.5j) implies 
that the deformed Laguerre polynomial lIC'\x, 7, 6) satisfies the three-term recurrence ()2.9()j) 
with 



an = 2n + a + l+ (n{2n + 1) + 2na + 7^5, (4.9) 
bn = n{n + a)(l + {2n + a- 1)6 + {{n - l){n + a) + 7)5') , c„ = (-l)"/n!, (4.10) 
(n + 1)4'J 1 (x, 7, 5) + - (2n + a + 1 + (n(2n + 1) + 2na + 7)^)) {x, 7, 6) 
+ {n + a)(^l + {2n + a- 1)6 + {{n -l){n + a)+ -f)6^^ L^^\{x, 7, 5) = . (4.11) 

In this case (j4.ip . the parameter ranges are 5, 7 > and a > — 1. So 6„ is positive. There- 
fore Llt\x,'y,6) is a set of orthogonal polynomials with respect to some positive definite 

inner product. From this three-term recurrence we obtain the differential equation for the 

00 

generating function G^"-* (t, x, 7, 6) = ^ t^L^n^ {x, 7, 6), 



d_ 



n=0 

a +- 



n-*t4 



t)G(")(t,x,7,5) = 0, (4.12) 



which is a third order linear differential equation with respect to t. 



Cr • For the solution {qj} of the Cr equation ()2.(i2|l . let us define 



Qj = y/ag^ Vj, 6 = —, a= — -l, (4.13) 

a 9s 

and introduce a degree r polynomial in x, having zeros at {yj}'- 

4-^^,5) = ^n(^-yi)- (4.14) 

This is a deformation of the Laguerre polynomial such that \imLl°'\x,6) = LI"\x), and 

5— >0 

obviously it is a special case of Li^\x,'y,6), 

4")(x,5) = 4")(x,0,5). (4.15) 
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The functional equation for Li^\x,6) is easily obtained from that of Li^\x,'-f,6) ()4.5j) 
and will not be presented. The three-term recurrence for {x, 6) reads 



[n 



+ 1)4+1 (a;, S)+ (x~{2n + a + l + n{2n + 2a + 1)5)) 4"^ {x, 6) 

+ (n + «) (l + (n - 1)5) (l + (n + a)6^L^^},{x, 6) = , (4.16) 



The value at the origin shows a simple deformation pattern 



4"H0,5)=(" + ")n(l+j5)- (4.17) 

Br : For the solution {qj} of the Br equation ()2.61|) . let us define 

Qj = VaglVj , S = —, a = — -1, (4.18) 

a 9l 

and introduce a degree r polynomial in x, having zeros at {Vj}'- 

4"H^,5) = ^ri(^-^i)- (4-19) 

This is a deformation of the Laguerre polynomial such that lim L^"-* (a;, 0) = LI^'^x), and 

<5— >0 

obviously it is a special case of Li^\x, 7, 6) 

Ll^\x,6)=Ll^\x,\{a + lf,6). (4.20) 

The functional equation of lI^\x,6) is easily obtained from that of 4"'' (a;, 7, 5) (j4.5p and 
will not be presented. The three-term recurrence for Zi"^ (x, 6) reads 

(n + l)Li'Ji (x, 5) + - (2n + a + 1 + (n(2n + 2a + 1) + i(« + 1)^)5)) 4") (x, 5) 

+ (n + «) (^1 + (n + ^izl)5)'L(f_\(x, (5) = 0. (4.21) 



Dr : As in the Calogero systems, the Dr is a special case (75 = of the Br theory described 
by Li ^\x, 6) = lI ^\x, 6), which has a zero at x = for all r. 
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4.2 Quadratic Confining Potential Case (2 or 3 Parameter Defor- 
mation) 

BCr- For the solution {qj} of the BCr equation ()2.72p . let us define 

/ X 90 ^ 91 + 92 ^ 9i92 . . 

Qj = VagoVj, d = —, e= -, a = 1, 7 = , (4.22) 

a a go 

and introduce a degree r polynomial in x, having zeros at {l/J}: 

L(")(x,7,5,e) = ^ri(x-yj). (4.23) 

It is a further deformation of the deformed Laguerre polynomial defined previously, 

lim L(°) (x, 7, 5, = 4°^ (a;, 7,5)- (4.24) 
The symmetry between the two parameters a 5 is expressed as 

4")(x,7,fc,e-^) = e'^L(")(e-V,7,5,e). (4.25) 

If we define Lr"^ (x, 7, *= f = 5, ^2 f = f ) *= (l+£:)''I^r"^ (y^^, 7, 5, e), then this symmetry 
is more manifest, Lr"^(x, 7, 5i, ^2) = -f'r°^(x, 7, ^2, ^i)- 

The equation for the equilibrium p.72|) can be written as (we replace r by n) 

" (^^. - r^5f -yl _ y, - + 1/. + y] + ^(« + l)v^ - 7^^ 

i-i(%+^v^)2-yfc' y,+i^ Vi-^Ts y=-^in 2/|-^(a + l)v^%-75' 

The equivalent functional equation for Li"'*(x, 7, 5, £:) in > 1) reads 

^ ( - + + (y' + + 1) v^y - 7'^) ((y + ^v^ 7, S, e) 

+ [y + ^^)[y-^^){y-^^)[y'-^ia + l)^rSy- 7^) L ^ ^(y - ^V~SY,^,S,e)'^ 

= 2{A^y' + B^y' + Cn)L^:\y^ 7, ^, (4-27) 

because the LHS is a degree 2n + 4 even polynomial in y with real coefficients which vanishes 
dX y = ±yj. Here An = Ai^\'y,6,e), Bn = Bli'\'y,S,e) and C„ = CiC\'y,S,e) are real 
constants: 

An = 1, (4.28) 
Bn = -{6-^ + 2n + a + - (2n + a + i) - (7 - ^^±i + 2n{n + a)) 5, (4.29) 

Cn = (a + 1 - 27 + (n + 7)5)^"' - ^{n + 7)5 - in(n + a)5^ (4.30) 
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This functional equation contains all the information of the equilibrium. In a proper limit 
£ — > 0, the above functional equation ()4.27|) reduces to that of Li"^(x,7,5) ()4.5p . This 
functional equation can be written as a difference equation as previous examples. 

The functional equation ()4.27|) implies (see Appendix) that the deformed Laguerre poly- 
nomial Ll^\x,'y,6,e) satisfies the three-term recurrence ()2.90|) with, 

o I 111 ^o + Xie + X2e'^ //I Qi\ 

ttn = 2n + a + l-\ -^—z , (4.31) 

"2n"2n+2 

bn = n{n + a){l + {2n + a - 1)5 + {{n -l){n + a)+ 7)^^) 

X (^1 + (1 + (n - l)^)^) (^1 + (2n + a - l)5e + {{n -l){n + a)+ 7)(5V) 

X;^^ , c„ = (-l)7r.!, (4.32) 

"2n-l"2n^2n+l 

dm = 1+ (l + {m + a-l)5^e. (4.33) 

Here Xq, Xi and X2 are 

Xq = [n{2n+l) + 2na + -i^5, (4.34) 

Xi = -(2n + a + 1) - (2n{n + a + 1) + (a + 1)^ - 27)5 

+ (^n(4n^ + (6a + l)n + 2a^ - 1) + {2n + a - 1)7) 5^, (4.35) 

X2 = -(2n + a + 1) - (gu"^ + 3{2a + l)n + 2a{a + 1) - 7)5 

-(^An^ + 3(2a + l)n^ + {Aa^ + 4a - l)n + (a - l)(a + 1)^ - (2n + a - 1)7)5^ 

+n(n + a) (^2n(n + a) - a - 1 + 27 j (5^ (4.36) 

Namely L^^^ {x, 7, 5, e) satisfies 

{n + l)4"^li(x, 7, 5, e) + (x - a„)4°)(x, 7, S,e) + - 6„ Lt\(x, 7, ^, e) = 0. (4.37) 

n 

In this case fl4.22|) . the parameter ranges are 5, e, 7 > and a > — 1. So 6„ is positive. From 

this three-term recurrence we obtain the differential equation for the generating function 

00 

n=0 

(2 ^ 2 

( ('^2n-l'i2n'^2n+l'^2n+2) „ (tt: + ^) ~ ('^2n-lC^2n'^2n+lC^2n+2fln) 

+ {- d2n-id2ld2n+id2n+2bn) t ) G^"^ (t, X, 7, (5, e) = , (4.38) 

which is an 8-th order linear differential equation with respect to t. The special case of 5 = 
gives the original Laguerre polynomial, l}n\x^ 7, 0, e) = (1 -|- £ )~"Ln"''((l -|- e )x). 
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Cr : For the solution {qj} of the Cr equation ()2.71|) . let us define 

Qj = y/agsyj, 5 = —, € = -, a= — -l, (4.39) 
a a gs 

and introduce a degree r polynomial in x, having zeros at {?/|} : 

i=i 

This is a further deformation of the deformed Laguerre polynomial defined previously such 
that lim4°)(a;,5,e) = 6) ()4.14|) . and obviously it is a special case of (x, 7, 5, e), 

e— >0 

Li''\x,5,e) = Li''\x,0,5,e). (4.41) 
The value at the origin shows a characteristic deformation pattern 

^" ^°'^''^=l n Jlll + (l + (a + ,+.)5).- ^'-"'^ 
Br : For the solution {qj} of the 5^ equation p.70|l . let us define 

Qj = VagZvj, 5 = —, e=-, « = — -!, (4.43) 
a a gi 

and introduce a degree r polynomial in x, having zeros at {?/|} : 

^l"H^,5,.)^=^^n(a:-y|). (4.44) 

This is a further deformation of the deformed Laguerre polynomial defined previously such 
that limL(°)(x,5,e) = L(°)(x, 6) ()4.19|1 . and obviously it is a special case of Lr (x, 7,5, e), 

(x, 5, e) =L(-)(x,i(a + 1)2,5,5). (4.45) 

Dr : As in the Calogero systems, the Dr is a special case (75 = of the Br theory described 
by Li ^\x, 6, e) = Lr ^\x, 6, e), which has a zero at x = for all r. 

Deformation of the identities : Before going to the systems with trigonometric poten- 
tials, let us present the one and two parameter deformation of the identities between the 
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Hermite and Laguerre polynomials, ()2.20|) and ()2.2H) . which could be considered as conse- 
quences of the folding ()2.5|) - ()2.7|) and ()2.50|) - ()2.52|) of the rational potentials. The one and 
two parameter deformation of the even degree identities ()2.2m) are 

2~^'H2r{x,6) = {-lYr\Li''^\x',6), (4.46) 
2-^^H2r{x,6,e) = {-iyr\Li'^'\x\6,e), (4.47) 

which are connected with the folding A2r-\ Cr ()2.5p . The one and two parameter defor- 
mation of the odd degree identities ()2.21|) are 

2-^''-'H2r+i{x,6) = x{-lYr\Li^\xM6), (4.48) 
2-^'-'H2r+i{x,6,e) = x{-lYr\Li^\xMS,e), (4.49) 

which are related to the folding 742,- — > BCr ()2.7|) . The one and two parameter deformation 
of the identities between Laguerre polynomials ()2.22|) are 

{r + l)L^Y+^{x,6) = -xU^\x,5), (4.50) 
(r + l)Li;l)(x,(5,£) = -xL«(x,5,£), (4.51) 

which are related to the folding Dj.+i Bj. (|2.6|) . 

5 Deformation of the Jacobi Polynomial 

5.1 Ar-l 

For the Ar-i case, the equilibrium position, ie the solution {qj} of the Ar-i equation ()2.83|) . 
is the same as the original Sutherland system ()2.31|) . Therefore the polynomial describing 
the equilibrium is same as the original Sutherland system, the Chebyshev polynomial of the 
first kind Tr{x) ()2.33p . In other words the Chebyshev polynomials are not deformed in the 
present scheme. 

5.2 B'Cr 

B'Cr : For the solution {qj} of the B'Cr equation ()2.88|) . let us define 

2 tanhgs tanh25(L t&nh2gL . 

= tanh gM, ot = 1 1, p = 1, (5.1) 

tanh gu 2 tanh gM 2 tanh gM 
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and introduce a degree r polynomial in x having zeros at {cos2gj} : 



P("'^)(x, 6) 2-^(°+^+20 n(:^ - cos2g,) . (5.2) 

It is a deformation of the Jacobi polynomial ()2.H4|1 such that 

limPj"'^)(x, 5) = Pj"'^)(a;)- (5.3) 

5— >0 

The equation for the equilibrium ()2.88|) can be written as (we replace r by n) 

irf cos 2qj + i sin 2qj — cos 2qk 

k=i cos 2qj - i sin 2qj - cos 2qk 
^ V6cos2qj + tsm2qj + ^ ("^ " /^)^ £l| + ^ 2(/3 + l)v/^ 

" v^cos2g,-^sin2g, ,i^-.(a-/5)v^ - ^2{f3 + 1)V6 ' 



Since qj can be restricted to < < 7r/2, sin2gj is sin2gj = yl — cos^ 2gj. From this 
equation, for — 1 < x < 1, we obtain the following functional equation for Pr"''^\x,6) 
{n>l), 



f3)V6il + x)+ zv^r^) (2{(3 + l)V6x + iVT 



x"^ 



+ x) — 



2iVl - x2 (A„x2 + P„x + C„)P^"''^)(x, 5), (5.5) 



because the LHS is — x"^ times a degree n + 2 polynomial in x with real coefficients 
which vanishes at x = cos2qj. Here A„ = Ai^'^\6), Bn = B^'^\6) and C„ = Cii"^\S) are 
real constants: 

A„ = -(l-5)-("-i) X i(^(i + (a-/5)v^)(l + 2(/3 + l)v^)(l + v^)2"-i 

+ (1 - (a - /?)v^)(l - 2{P + l)v^)(l - v^)2"-i), (5.6) 
P„ = -(a-/3)(l + 2/?)5, (5.7) 
C„ = (l-5)-"x i((l + (a-/3)v^)(l + 2(/3+l)v^)(l + v^)^"-^ 
+(1 - (a - P)V6){1 - 2(/3 + l)v^)(l - v^)2"-^ 
-2(a-/?- 1)(1 + 2/3)5 (1-5)"-^). (5.8) 
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The functional equation ()5.5|) contains all the information of the equilibrium. In the 6—^0 
limit, this functional equation reduces to the differential equation of the Jacobi polynomial 

(EMD. 

The three-term recurrence fl2.90|) for the deformed Jacobi polynomial Pli''^\x,6) is a 
consequence of the functional equation ()5.5|) (see Appendix). The constants in ()2.90|) are 

an = 

xi((l + 2(3) ((1 + (a - f3)V~5){l + 2(/3 + l)v^)(l + v^)^"-^ 
+(1 - (a - /3)v^)(l - 2(/5 + l)v^)(l - v^)2"-i) 
+2(a-^-l)(l-4(/3+l)25)(l-5)"-i) X — ^ , (5.9) 



^((i + v^r-(i-v^r 



d2n C^2n+2 



2^5 



:^((1 + 2(/3 + l)v^)(l + v^)"-i + (1 - 2(/5 + l)v^)(l - v^)""^) 
;^((1 + 2(/? + l)v^)(l + v^)2"-2 - (1 - 2{(3 + l)v^)(l - v^)'^ 
:i((l + (a - /3)v^)(l + v^)"-i + (!_(«- (3)V~5){1 - 75)"-^) 



2n-2 



^((1 + (a - /?)v^)^(l + 2(/? + l)v^)(l + v^) 

-(1 -{a- /3)v^)'(l - 2(/? + l)v^)(l - v^)2"-2 
3 ^^7^' Cn = 2-\-^^:'-), (5.10) 

"2n-l "2n "2 



''2n "2n+l 

= ^('(l + (a-/5)v^)(l + 2(/3 + l)v^)(l + v^)™-2 
2v 5 ^ 



-(1 - (a - /?)v^)(l - 2(/? + l)v^)(l - 75)™-'). (5.11) 
Namely P^'^^ (x, 5) satisfies 



2(n + l)(a + /3 + n + l) 



-Pi:f (x,5) + (x-a„)P^^)(x,5) 



{a + 13 + 2n + l){a + 13 + 2n + 2) 

^^''"%„Prf (x,5) = 0. (5.12) 



(a + /3 + 2n- l)(a + /3 + 2n) ^ 



2n{a + 13 + n) 

In this case (15.11) . the parameter ranges are 5 > and a > /? > — 1. So 6„ is positive. From 
this three-term recurrence we obtain the difference equation for the generating function 
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n=0 



f 2 ^ 2 

( (c?2n-l'^2n'^2n+l'^2n+2) ^ (^TT + x) — ((i2n-l'^2n'^2n+lC^2n+2fl 



n^iA dt 



+ {d2n-ld,ld2n+ld2n+2bn) , t ) G^i (t, X, 5) =0. (5.13) 

at 



It is interesting to note that the three-term recurrence ()5.9|1 - ()5.10|1 simphfies drastically 
in a 'strong' coupling limit, 



qm +00 5-^1, ql, gs ■ fixed, 



ao 



(3 — a 



a) (2/5 + 1) 



a + p + (a + /? + 2)(a-/3 + l)' 



a„ = 0, {n> 2), 
a + /3 + 2 



, _ 4(/3 + l)(a + l + (/3 + l)(a-/3)^) ^ ^ 

' (a-/3+l)(a + /5 + 2)2(2/3 + 3) ' ' 2{a - p + l){2p + 3)' 

Cr : For the solution {qj} of the Cr equation (j2.86p . let us define 



(5.14) 
(5.15) 

(n > 3). (5. 16) 



6 = tanh^ gs, a 



tanh 2gL 



- 1, 



2 tanh gs 

and introduce a degree r polynomial in x having zeros at {cos2gj} : 

6) 2^'("-J+0 f[{x - cos2g,) . 

It is a deformation of the Gegenbauer polynomial ()2.35|) such that 

limd"^^^(a;,5)=d""'^\x), 

<5^0 



(5.17) 



(5.18) 



(5.19) 



and obviously it is a special case of Pr"'^\x, 6) with definite parity, 

d"^^\x,5) = ~V,("'")(x, 5), d"^^Va^,<^) = (-l)'Cr"^^\a;,5). (5.20) 



The functional equation of Cr^^^'{x,S) is easily obtained from that of Pr^'^\x,6) and will 
not be presented here. The three-term recurrence for clC~^^\x,6) reads 

^Ji±i-^c'::,%,s)+.c:r-%,s) ^?^^±^i„ct-%,s)=o. (5.21) 



n 
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B'j. : For the solution {qj} of the B'^ equation ()2.85|) . let us define 



, 9 tanh oq 

S = tanhV, a = - 1, 5.22 

tanh ql 

and introduce a degree r polynomial in x having zeros at {cos2gj} : 

r 

P,(-)(x, 6) 2-("-;.+2^) - cos2g,) . (5.23) 
Obviously it is a special case 

Pj")(x,5) =Pj"'-^)(a;,5)- (5.24) 

The functional equation and three-term recurrence are obtained from those of 

P^^H^,^)- 

5.3 

CI : For the solution {qj} of the equation ()2.87|) . let us define 

2 tanhgf^, 

= tanh gs, a = 1, (5.25) 

tanh gs 

and introduce a degree r polynomial in x having zeros at {cos2gj}: 

r 

Ci"+^)(x,5) = 2'-("-J+'-) l[{x - cos2g,) . (5.26) 

It is a deformation of the Gegenbauer polynomial ()2.35|) with definite parity 

\imCi'"^'^\x,6) =d'"^'^\x), ci'"^'^\-x,6) = {-lYcl'"^'^\x,6). (5.27) 

<5— >0 

The functional equation for C^ri^~^'^\x,5) {n > 1) reads 

(VSx- iVr^^) ((a + 1) v^x + iVl^Jc^"^'^'^ {t^6^ + ^f^v^T^^, 5 
-(ySx + zv^r^^) ((a + l)v^a; - tVl^J C^"'^'^^ (irf ^ - zf^TT^^, 5 
= 2iVl - x2 (A„x2 + P„x + C^)^™"^^^^, 5). (5.2^ 

Here = A^n^^'\6), P„ = Pi"^^\(5) and C„ = Cn^^'\6) are real constants: 



A„ = -(l-5)-("-i)xi((l + (a + l)v^)2(l + V5)2-i 
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+ (1 - (a + l)v^)'(l - 75)'"-^], (5.29) 
Bn = 0, (5.30) 
Cn = (l-5)-"x i((l + (a + l)v^)2(l + v^)2-i 

+ (1 - (a + l)v^)'(l - 75)'"-^ - 2a^6 (1 - 5)"-^). (5.31) 

This functional equation contains all the information of the equilibrium. 

The deformed Gegenbauer polynomial C^^^\x,6) satisfies the three-term recurrence 
fimn|l (see Appendix) with 

an = 0, (5.32) 



"2^ 



^^((1 + {a + 1)V6){1 + V6r~' + (!-(« + l)v^)(l - v^)"-i 



■ . . . ^n = 2"rM, (5.33) 

"2n-l "2n+l 

dm = -^({l + ia + l)V5f{l + V5r'^-{l-{a + l)V5Yil-Vsr-A. (5.34) 
2v5 ^ ^ 

Namely ci^^^^x, 5) satisfies 

" + ' -Ci-^^ (X, S) + X C'r-' (X, 6) + ^- + ^—' (-, ^) = 0. (5.35) 



2a + 2n + l ^ ' ^ ^ ' ^ n 

In this case ()5.25p . the parameters are 6 > and a > —1. So 6„ is positive. From 
this three-term recurrence we obtain the difference equation for the generating function 

oo 

G("+^)(t, X, 6) = ^rC'i""*"^V, in a similar way to that of Pt''^\x, 6). 

n=0 

In a 'strong' coupling limit 

gs +00 5 I, Ql- fixed, (5.36) 
the three-term recurrence (|5.33p - (|5.34p simplifies drastically: 

5.4 5^ 

: For the solution {qj} of the equation ()2.84|) . let us define 

„ 2tanh2f 

5 = tanh^c/i, a = — ^ - 1, 5.38 

tanh^ii 
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and introduce a degree r polynomial in x having zeros at {cos2gj}: 

r 

Pl:^\x,5) = 2-\-~\+^'') - cos2g,) . (5.39) 

It is a deformation of the Jacobi polynomial ()2.34|1 such that 

limPj")(x,(5) = Pj"'-^)(x). (5.40) 
The functional equation for P^\x,5) [n > 1) reads 

(VSx- iVT^'^ (^v^ (1 + x) + iVT^JpjC^ (t5 a; + ij^VT 



- (y6x + iVl^^ (1 + x) - iVl^^ P^") {tS^- «fr|v^l^^, S 

= 2iVl-x^ {Anx^ + Bnx + Cn)PlC\x,5). (5.41) 

Here = aI^\6), Bn = B^\6) and C„ = C^\6) are real constants: 

A„ = -(l-5)-(«-l) X i((l + l(« + 1)75)2(1 + V^)2«-l 

+ (1 - i(« + l)v^)'(l - 75)'"-^), (5.42) 
Bn = i(l-a2)5, (5.43) 
C„ = (l-5)-"xi((l + l(a + l)v^)2(l + v^)2""i 

+(l-i(a + l)v^)2(l- v^)2"-i 

-(1 + _ 1 (tt + 1)2^)5 (1 _ s )"-i) . (5.44) 

This functional equation contains all the information of the equilibrium. 

The deformed Jacobi polynomial pl"^ (x, 6) satisfies the three-term recurrence ()2.90|) (see 
Appendix) with 



an = il-a')il-6r'4^^^, (5.45) 



^2n W2n.+2 



bn = 4-^((l + v^)"-^-(l-v^)"-^) x-i=((l + v^)"-(l-v^)") 



2v^V' ' y 2y/S 

rl'4: rl rl 
d2n-l C^2n 'J^2n+1 



2-"r„ ), (5.46) 



c^n^ = ^((l + i(a + l)v^)'(l + v^)™-'-(l-|(a + l)v^)'(l-v^)™-'), (5.47) 



dm 



i((l + i(a + l)v^)(l + V~5r-^ + (1 - i(a + l)v^)(l - v^)'"-^). (5.48) 
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Namely Pn"\x,6) satisfies 



2{n + l)ia + n) ^(.) . . .p^a). ^^ 

(a + 2n-2)(a + 2r2- 1) ^ 



. , ^ bnPr-\ix,6) = 0, (5.49) 

2ri(a + n — 1) 

In this case fl5.38p . the parameters are 6 > and a > — 1. So 6„ is positive. From this three- 
term recurrence we obtain the difference equation for the generating function G^"\t, x, 6) = 

oo 

J2 t'^Pn\^^ ^) in a similar way to that of Pt'^\x, 6). 

n=0 

In a 'strong' coupling limit 

gL +00 5 1, gs ■ fixed, (5.50) 
the three-term recurrence fj5.45|) - ()5.48|) simplifies drastically: 

Oo = -1, ai = a„ = 0, (n > 2), (5.51) 

^1 = 0' ^' = T^^^ ^" = 7' (^^3). (5.52) 
(a + 6)'^ 4 

: As in the Sutherland systems, the Dr is a special case gs = of the Bj. theory 
described by Pr~^\x, 5), which has a zero at x = ±1 for r > 2. 

Deformation of the identities : Before closing this section let us briefiy discuss the 
deformation of the identities between the Chebyshev and Jacobi polynomials ()2. 37^ - 1)2. 38p 
and those between the Jacobi polynomials ()2.39|) . As remarked in section I^TT| the Chebyshev 
polynomials describing the equilibrium of the Ar-i systems are not deformed in our scheme. 
Therefore no deformation of the identities ()2.37|) - ()2.38|) exists. The folding Dr+i Br ()2.52p 
leads to the identity between the deformed Jacobi polynomials Pr""* {x, 6) ()5.39|) associated 
with the Br systems 

2(r + l)P,y^(a;, S) = r{x - l)P!:^\x, S) , (5.53) 

which is a deformation of the identity ()2.39p . As remarked at the end of section I^Tl we 
have not been able to deform the identities between the Gegenbauer and Jacobi polynomials 
()2.40|) - ()2.41|) as they do not seem to have a root theoretic explanation. 

Among the deformed Jacobi polynomials P!f''^\x,6) for various a and (3, two cases (i) 
a = f3 = —1/2, (ii) a = —(3 =1/2 are not deformed: 

p(-V2,-i/2)(^^5^ = Pt'''-"'\xl Pi'''-'''\x,6) = Pi'/''-'/'\x). (5.54) 
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The first is the Chebyshev polynomial of the first kind Tn{x) oc cos rvf, x = cos (p, as remarked 

('1/2 1/2) 

in section I5.ll The second case is Pn (x) oc sin((2n + l)(p/2)/ sm((p/2), x = coscp. 

In both cases, the zeros of Pn"''^\x) are equally spaced. There is a third case ^ of 
equally spaced zeros of the Jacobi polynomial Pn"'^\x), for a = /5 = 1/2, corresponding to 
the Chebyshev polynomial of the second kind Un{x) oc sinny?/ simp, x = cosy), which is, 
interestingly, deformed. We have no explanation to offer. 

6 Summary and Comments 

We have derived certain deformation of the classical orthogonal polynomials (the Hermite, 
Laguerre, Gegenbauer and Jacobi) describing the equilibrium of a class of multiparticle 
dynamics, the Ruijsenaars- Schneider systems. The R-S systems are 'good' deformation of 
the Calogero and Sutherland systems whose equilibrium points are described by the zeros 
of the above classical orthogonal polynomials. As remarked in the text these deformed 
polynomials do not belong to the g-deformed orthogonal polynomials ^3] or their analogs 



The quality and quantity of the knowledge of these new polynomials are rather varied. 
The one parameter deformation of the Hermite polynomials section ITT] is best understood. 
Its three-term recurrence (13.111) tells that it is the simplest possible deformation which re- 
duces to the original Hermite polynomial without rescaling etc in the zero deformation limit 
{6 —>■ 0). As shown in some detail in section ITTl the generating function ()3.14|) and the 
weight function ()3.21|) are known explicitly. Some identities connecting the Hermite and 
Laguerre polynomials ^TKi^ - ^nTli are nicely deformed (I0H|) - (I07I) . (IOHl) - (IOni) . It is 
interesting to note that some non-trivial identities of the Hermite polynomials are preserved 
after deformation. For example, the 'addition theorem' reads 

E ^^H^,{x,)---H^^{xJ=^-^hJ--), (6.1) 

rii'. ■ ■ ■ rim'. nl V \a\ J 

TilH [-nm—n 

in which the following notation is used: a = *(ai, ■ ■ ■ , am), x = ^{xi, ■ ■ ■ ,Xm), |a| = 



>y + ■ ■ ■ + a^, a ■ X = aiXi + ■ ■ ■ + CKrnXm- The deformed version is 



00 



<^H„,(.,i,)...H„.,.(x„ A) = H:;^„(1^ * ). (6,2) 
, Till- ■■rim' ^ oti/ V ai,/ n\ V \a\ \a'^ 

*^l,'")"m — III 
n-i-\ \-nm — 'n 
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Both can be derived from the generating functions. 

After a long and laborious search through existing literature, we find out that all the 
deformed orthogonal polynomials introduced in the main text can be related to particular 
members of the Askey-scheme of hypergeometric orthogonal polynomials 16 . 

The deformed Hermite polynomial Hn{x,6) is related to the Meixner-Pollaczek polyno- 
mial pk^\x-(t)) (§1.7 of nni): 

i/„(x,5)=n!v^"pi^^(^;|). (6.3) 

This identification allows a simple expression of its general term in terms of a (truncated) 
hypergeometric series 2F1: 

-.n/6\ „/-n,j + i3f 



6 

^«"g,-if(^)n(..-L,,i,).n(-L,f, 



^- (6,4) 



k=0 ^ ' j=0 ^ " j=k 

which is deformation of ()2.18j] . The two-parameter deformation of the Hermite polynomial 
Hn{x,6,e) ()3.27|) is a special case of the continuous Hahn polynomial (§1.4 of [T^): 

''-^"'^'^^= (n-l + | + ^)/ 47!'^'fe-^-fej- (^-^^ 

The (two-parameter) deformed Laguerre polynomial lI^\x, 7, 6) ()4.2|) is the continuous dual 
Hahn polynomial (§1.3 of [16j) with rescaling: 

7, 5) = (y 5 ^ ' "1 + 1' "2 + 1) , (6-6) 

in which ai and 02 are the two roots of — (a — l)x -|- 7 — a = 0. The (three-parameter) 
deformed Laguerre polynomial Li^\x,'j,6,e) ()4.23j) is the Wilson polynomial (§1.1 of [TT)] ) 
with rescaling: 

5" fy^ 1 1 

The deformed Jacobi polynomial Pn"''^\x,6) ()5.2|) is a special case of the Askey- Wilson 
polynomial (§3.1 of [TB]): 

P,!"'''>(rc,^) = 2-'^"C>+'J«")(a62,»-l;,);;'p„(i;a,6,-6,-l|9), (6.8) 

1 + 75 l + (a-/3)\/5 l + 2(/3 + l)v^ 
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4")(y^ 7, 5, e) = ^ M ; t , ;r ' «i + 1' «2 + 1) • (6-'^) 

n! (n + a + i + TT)n V 5 d / 



The deformed Gegenbauer polynomial ci^^^^x, 6) ()5.26p is again a special case of the Askey- 
Wilson polynomial (§3.1 of [IF]): 

Ci''^'^\x,6) = ("i+")(aV-';g);Vn(a;;a,a,-a,-a|g), (6.10) 

g=i^ = e--, a^ = i^i^±M = e--. (6.11) 
1 + V5 1 + (a + 1)75 

Another deformation of the Jacobi polynomial P-jf^ (x, 6) (j5.39p is also a special case of the 
Askey- Wilson polynomial (§3.1 of pT)]): 

l-l(a + l)v/^ 

g = ^ = e a = 7 ^ = e (6.13) 

1 + 75 l + i(a + l)v^ ^ ' 



In all these formulas the Pochhammer symbol {a)k = nJ=o('^+i) g-extension (a; g) 



nJ=o(l ~ (^1'') are used. 

Various 'strong' coupling limits ()5.14|) . ()5.36|) and ()5.50|) of the deformed Jacobi type 
polynomials Pn°''^\x,6), cli'~^^\x,6) and P^\x,6) simply correspond to the ^crystaV limit 
g ^ 0+ of the Askey- Wilson polynomials ^H^i and (jfTT^ . 

For all these polynomials discussed in the present paper, the general term can be expressed 
in terms of various hypergeometric functions 2-^1, 3-^2, and 4^3. A Rodrigue type formula, 
the generating function and the weight function etc can be written down by using the general 
formulas of the Askey-scheme of hypergeometric orthogonal polynomials [T^ . 

As remarked repeatedly in the text, the equations determining the equilibrium positions, 
(Pl)n|l - (ITM|1 . ^IW^ - ^rn^ and drrpi - ipipi look similar to the Bethe ansatz equation. For 
the simplest spin 1/2 XXX chain with sites, the Bethe ansatz equation reads 

ri^^^^=(2i±^)", (.=1,...,/), (6.14) 

in which / < is the number of up (down) spins. This looks similar to the rational A 
type equations ()2.60|) and ()2.69|) with a special choice of the potential w{x) function, the 
A^-th power rather than linear or quadratic. The corresponding functional relation, Baxter's 
t-Q relation, looks very much like the functional equations (13. 6j) and (j3.31|) . It would be 
interesting to pursue the analogy further. 
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Appendix: Relation between the Functional Equation 
and the Three Term Recurrence 

In this appendix we show the relation between the functional equation and the three-term 
recurrence, without proof. Since the normalisation of polynomials is irrelevant, we consider 
monic polynomials fn{x) with real coefficients and without superscript 'monic'. 
The three-term recurrence of fn{x) ()2.90p is 

f n+lix) = {X - an) f nix) -bnfn~lix) , (u > 0) , (A.l) 

with f-i{x) = and fo{x) = 1. The explicit forms of an and 6„ can be read from ()3.1()|1 . 

(uni, ra-iniD, ^M-^M, ra-dnni), ^m-^m, ^m-^m- 

The functional equations for the deformed Hermite, Laguerre and Jacobi polynomials 
have the following forms {n > —1): 

Hermite : h{x) fn{x + iVd ) + e h{x)* fn{x ~ iVd ) = 2 i~ gn{x) fn{x) , (A. 2) 

Laguerre : h{y)fn[{y + ^V6)')+eh{y)*u(^{y-^^/6)^') = 2 g^iy) Uv') , (A.3) 
Jacobi : /^(x)/„(i±fx + zf^v/r^) +e/i(x)7„(i±fx-2f^v/r^) 

= 2(zv/r^) ' gn{x)fn{x), (A.4) 

where e is 

e = -l : i^„(x,5),4")(x,7,5),4"Ha;,5),Li")(x,5),P^^)(x,5),C7i"+^\x,5)^ 
6 = 1 : Hnix,S,e),Ll:\x,^,S,e),Ll^"\x,5,e),Ll:\x,5,e),c!C^'^\x,S), p(") {x,S). 

The explicit forms of h{x) and gn{x) can be read from (ICTIl . (H3|) . KT7\ . ^^^M 

and (ICTll ^. 

As the first step, we show the following property of the functional equations: 
^The explicit forms of the function h{x) are derived from the equations for the equihbrium. For the 
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Proposition A.l The solution of the functional equation ()A.2|) - ()A.4|) . if exists, is unique 
up to an overall normalisation. 



By this proposition, it is sufficient to construct one solution of each functional equation 
explicitly. We will do this by using the three-term recurrence. 
As the second step, we show the following relation: 



Proposition A. 2 Three-term recurrence flA.l|) implies the relation {n > 0)^ 

Hermite : i^V6 h{x)fn{x + iV6) 

{x - a„)(5(„+i(x) - gnix)) + iVd 9n{x) ) fn{x) 

-bn(^9n+l{x) - gn-l{x)^ fn-l{x) , (A. 5) 



Laguerre : 2i 2 V6 yh{y) fn(^{y + iVS Y 

= ({y^ - an){gn+i{y) - gn{y)) + 5gn{y) + "^iVs ygn{y)^ fn{y^) 

-bnign+iiy) - gn-iiy))fn~iiy'^) , (A.6) 
. . 1+1 . . 

Jacobi: UVT^) ' M/i(x)/„fi±|x + zf^yr^^) 

(x - an)gn+i{x) - ij^x- an)gn{x) + 2iVl - x"^ ^ gn{x)^ fn{x) 



-hn{gn+l{x) - gn-l{x))fn~l{x) ■ (A. 7) 

In the proof of this proposition by induction, we encounter the equation, 

X„(x)/„(x)-F„(x)6„/„_i(x)=0 or Xn{y)fn{y^)~Yn{y)hnfn-i{y'') = ^- (A.8) 

Here Xn{x) and Yn{x) are as follows: 
Hermite: 

Xn{x) = (x - anf{gn+i{x) - gn[x)) - (x - a„)(x - a„+i)(5f„+2(x) - gn+\{x)) - 6gn{x) 

-hn{gn{x) - gn^2{x)) + &n+l (5'n+2 (a^) " gn{x)) , (A. 9) 
K„(x) = {X - an){gn+l{x) - gn~l{x)) - {X - an-l){gn~l{x) - gn-2{x)) 

-(X - an+l){gn+2{x) - gn+l{x)), (A.IO) 



deformed Hermite polynomials, gn{x) are determined by the consistency of this functional equation. For the 
deformed Laguerre and Jacobi polynomials, however, gn{x) are not determined uniquely by these functional 
equations. We have fixed gnix) by using some empirical knowledge of their lower degree members. 
•^Although a„ vanishes for the deformed Hermite cases, we keep a„ in this generic formula. 
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Laguerre : 

Xn{y) = {y^ - an - S) ( (y^ - a„)(5'„+i(y) - gn{y)) + Sgn{y) 



-bn (gn iy) - gn-2{y)) + hn+i ign+2 {y) - gniy)) , ( A. 1 1 ) 

Yn{.y) = iy"^ - an - 26){gn+iiy) - gn-i{y)) - iy'^ - an-i){gn-iiy) - gn-2{y)) 

-{y^ - an+i){gn+2{y) - gn+i{y)) , (A. 12) 

Jacobi : 

Xn{x) = ijZ^X - an)({x - an)gn+l{x) - ij^X - an)gn{x) 



-{x- an) (^{x - an+i)g„+2{x) - X - an+i)gn+i{x)^ - (f^)^(l - x^)gn{x) 

-K{gn{x) - gn~2{x)) + hn+l{gn+2{x) - gn{x)) , (A.13) 

Yn{,x) = (Yrfa; - an)(5'n+i(a;) - gn-i{x)) + (x - an-i)gn-2ix) - ij^x- an-i)gn-iix) 

-(x - an+i)gn+2{x) + X - an+i)gn+i{x) . (A.14) 

It is easy to see Xn{x) = Yn{x) = by using the explicit forms of a„, 6„ and gn{x). 
As the third step, we show the following: 

Proposition A. 3 The polynomial defined by the three-term recurrence ()A.1|) satisfies the 
functional equation ()A.2|) - ()A.4|) . 



Therefore we obtain the following: 

Proposition A. 4 The polynomial defined by the functional equation ()A.2|) - ()A.4|) satisfies 
the three-term recurrence ()A.1|) . 

References 

[1] F. Calogero, "Solution of the one-dimensional A^-body problem with quadratic and/or 
inversely quadratic pair potentials", J. Math. Phys. 12 (1971) 419-436. 

[2] B. Sutherland, "Exact results for a quantum many-body problem in one- dimension. II", 
Phys. Rev. A5 (1972) 1372-1376. 



43 



[3] J. Moser, "Three integrable Hamiltonian systems connected with isospectral deforma- 
tions", Adv. Math. 16 (1975) 197-220; J. Moser, "Integrable systems of non-hnear 
evolution equations", in Dynamical Systems, Theory and Applications] J. Moser, ed.. 
Lecture Notes in Physics 38 (1975), Springer- Verlag; F. Calogero, C. Marchioro and 
O. Ragnisco, "Exact solution of the classical and quantal one-dimensional many body 
problems with the two body potential Va{x) = g^a^ j sinh^ ax", Lett. Nuovo Cim. 13 
(1975) 383-387; F. Calogero, "Exactly solvable one-dimensional many body problems", 
Lett. Nuovo Cim. 13 (1975) 411-416. 

[4] E. Corrigan and R. Sasaki, "Quantum vs Classical Integrability in Calogero-Moser Sys- 
tems", J. Phys. A35 (2002) 7017-7062. 

[5] O. Ragnisco and R. Sasaki, "Quantum vs Classical Integrability in Ruijsenaars- 
Schneider Systems," J. Phys. A37 (2004) 469-479; R. Sasaki, "Quantum vs Classi- 
cal Integrability in Ruijsenaars-Schneider and Calogero-Moser Systems", Czechoslovak 
Journal of Physics, 53 (2003) 1111-1118. 

[6] S. N. M Ruijsenaars and H. Schneider, "A New Class Of Integrable Systems And Its Re- 
lation To Solitons," Annals Phys. 170 (1986) 370-405; S. N. M Ruijsenaars, "Complete 
Integrability of Relativistic Calogero-Moser Systems And Elliptic Function Identities," 
Comm. Math. Phys. 110 (1987) 191-213. 

[7] J. F. van Diejen, "The relativistic Calogero model in an external field," 
Isolv-int/ 9509002 ; "Multivariable continuous Hahn and Wilson polynomials related 
to integrable difference systems", J. Phys. A28 (1995) L369-L374. 

[8] J. F. van Diejen, "Difference Calogero-Moser systems and finite Toda chains" , J. Math. 
Phys. 36 (1995) 1299-1323. 

[9] J. F. van Diejen, "Integrability of difference Calogero-Moser systems" , J. Math. Phys. 
35 (1994) 2983-3004. 

[10] F. Calogero, "On the zeros of the classical polynomials", Lett. Nuovo Cim. 19 (1977) 
505-507; "Equilibrium configuration of one-dimensional many-body problems with 
quadratic and inverse quadratic pair potentials", Lett. Nuovo Cim. 22 (1977) 251-253. 



44 



[11] S. Odake and R. Sasaki, "Polynomials Associated with Equilibrium Positions in 
Calogero-Moser Systems," J. Phys. A35 (2002) 8283-8314. 

[12] H. Awata, Y. Matsuo, S. Odake and J. Shiraishi, "Excited states of Calogero- Sutherland 
model and singular vectors of the Wn algebra," Nucl. Phys. B449 (1995) 347-374. 

[13] R.P. Stanley, "Some Combinatorial Properties of Jack Symmetric Functions", Adv. 
Math. 77 (1989) 76-115. 

[14] J. Shiraishi, H. Kubo, H. Awata and S. Odake, "A Quantum Deformation of the Virasoro 
Algebra and the Macdonald Symmetric Functions", Lett. Math. Phys. 38 (1996) 33-51; 
H. Awata, H. Kubo, S. Odake and J. Shiraishi, "Quantum Wn Algebras and Macdonald 
Polynomials", Comm. Math. Phys. 179 (1996) 401-416. 

[15] G. E. Andrews, R. Askey and R. Roy, Special Functions''^ , Encyclopedia of mathematics 
and its applications, Cambridge, (1999). 

[16] R. Koekoek and R. F. Swartouw, "The Askey-scheme of hypergeometric orthogonal 
polynomials and its g-analogue" , |math.CA/9602214l 

[17] A. J. Bordner, N. S. Manton and R. Sasaki, "Calogero-Moser models V: Supersymmetry 
and Quantum Lax Pair", Prog. Theor. Phys. 103 (2000) 463-487; S. P. Khastgir, A. 
J. Pocklington and R. Sasaki, "Quantum Calogero-Moser Models: Integrability for all 
Root Systems", J. Phys. A33 (2000) 9033-9064. 

[18] G. Szego, Orthogonal polynomials''^ Amer. Math. Soc. New York (1939). 

[19] 1. G. Macdonald, ^''Symmetric Functions and Hall Polynomials" , 2nd ed., Oxford Uni- 
versity Press, (1995). 

[20] T. S. Chihara, ^^An Introduction to Orthogonal Polynomials" , Gordon and Breach, New 
York, (1978). 

[21] S. Jing and W. Yang, 'A new kind of deformed Hermite polynomials and its application', 
|niath-ph/02 12011] 



45 



